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COSMIC RAY CURRENT-DRIVEN TURBULENCE AND MEAN-FIELD DYNAMO EFFECT 
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ABSTRACT 

We show that an a effect is driven by the cosmic ray Bell instability exciting left-right asymmetric turbu- 
lence. Alfven waves of a preferred polarization have maximally helical motion, because the transverse motion 
of each mode is parallel to its curl. We show how large-scale Alfven modes, when rendered unstable by cosmic 
ray streaming, can create new net flux over any finite region, in the direction of the original large-scale field. 
We perform direct numerical simulations (DNS) of an MHD fluid with a forced cosmic ray current and use the 
test-field method to determine the a effect and the turbulent magnetic diffusivity. As follows from DNS, the 
dynamics of the instability has the following stages: (i) in the early stage, the small-scale Bell instability that 
results in a production of small-scale turbulence is excited; (ii) in the intermediate stage, there is formation of 
larger-scale magnetic structures; (iii) finally, quasi-stationary large-scale turbulence is formed at a growth rate 
that is comparable to that expected from the dynamo instability, but its amplitude over much longer timescales 
remains unclear. The results of DNS are in good agreement with the theoretical estimates. 

It is suggested that this dynamo is what gives weakly magnetized relativistic shocks such as those from 
gamma ray bursts a macroscopic correlation length. It may also be important for large-scale magnetic field 
amplification associated with cosmic ray production and diffusive shock acceleration in supernova remnants 
(SNR) and blast waves from gamma ray bursts. Magnetic field amplification by Bell turbulence in SNR is 
found to be significant, but it is limited owing to the finite time available to the super-Alfvenicly expanding 
remnant. The effectiveness of the mechanisms is shown to be dependent on the shock velocity. Limits on 
magnetic field growth in longer-lived systems, such as the Galaxy and unconfined intergalactic cosmic rays is 
also discussed. 

Subject headings: ISM: cosmic rays - instabilities - magnetic fields - turbulence 



1. INTRODUCTION 

Astrophysical blast waves are strongly suspected of am- 
plifying the ambient magnetic field into which they propa- 
gate. Supernova remnants (SNR), given detailed models for 
their ultra high energy gamma ray emission, indicate mag- 
netic fields that are considerably stronger than the several /iG 
fields that are present in the interstellar medium. The exact 
strength of SNR magnetic fields depends on how small-scale 
bright spots are interpreted (Pohl 2009). 

Gamma ray bursts (GRB) afterglows, which are attributed 
to relativistic blast waves, are generally best fit with a mag- 
netic field strength that is much higher than interstellar mag- 
netic fields. It has been suggested that the Weibel instability 
is responsible for the magnetic field production/amplification 
(Medvedev & Loeb 1999), but several difficulties with this 
proposal (Blandford & Eichler 1987) remain unsolved. First, 
the fastest growing Weibel-unstable modes are of very small 
scale, the ion plasma skin depth, and they should decay away 
resistively over the hydrodynamical timescale of the blast 
wave. Second, the electrons in the actual interstellar medium 
are already magnetized as they enter the shock, and they 
should therefore freeze the magnetic flux. The Weibel insta- 
bility, which creates new flux, should therefore be suppressed 
at finite amplitudes, despite being unstable at infinitesimal 
amplitudes. It is not clear that these problems are resolved 
by simulations, which cannot be run over hydrodynamical 
timescales, and which do not always include initial electron 
magnetization. In this regard, the Bell instability (Bell 2004), 



which treats the thermal plasma as a magnetized fluid may be 
more relevant than kinetic approaches that ignore the electron 
magnetization, i.e., that their gyro radii are small compared 
with the relevant length scales of the system. 

Cosmic ray (CR) protons above 200 TeV and iron nuclei 
above 3 PeV are difficult to account for with standard SNR 
parameters (Lagage & Cesarsky 1983), and magnetic field 
amplification would solve this problem if it occurs on a large 
scale. The smaller the scale of field amplification, the lower 
the maximum energy of the CRs that can be accelerated by 
the shock (Eichler & Pohl 201 1). 

Simulations of magnetic fields in the presence of cosmic 
rays (Lucek & Bell 2000; Bell 2004) show magnetic field 
stretching. This is accompanied by a jumbling of the field 
lines into a more complicated geometry, and a smaller coher- 
ence length. It has remained unclear whether this is merely 
turbulent field line stretching or whether there is an additional 
mechanism. 

It has recently been noted by Bykov et al. (201 la) that a cir- 
cularly polarized Alfven wave gives a net electromotive force 
(EMF) along the direction of the original magnetic field. Be- 
cause cosmic ray protons preferentially excite Alfven waves 
of a particular circular polarization, they generate a net EMF 
along the original magnetic field. The field growth is given 
by the curl of this EMF which, for a plane wave, is kxB, 
yielding a growing Alfven mode, whose polarization vector 
k x B is perpendicular to the original magnetic field. To low- 
est order, this does not amplify the field but merely bends it, 
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and the question still remains as to whether (i) the field lines 
are merely getting stretched on this large scale, which would 
leave the net flux through any large-scale surface unchanged, 
or whether (ii) there is organized amplification, whereby the 
flux of magnetic field lines though a given large-scale surface 
is increased. 

In this paper we note that the situation can result in an a 2 
dynamo, in which a large-scale magnetic field can grow along 
its original direction and thus get amplified. The effect results 
from a nonlinear coupling between the individual waves that 
were demonstrated by Bykov et al. (201 la) to grow. However, 
there is a maximum size over which such amplification can 
be effective. As the quantity a has dimensions of velocity (in 
contrast to the stretching timescale, which has dimensions of 
time and which can be as short as the eddy turnover time of the 
turbulence), this maximum spatial scale is of the order of a T, 
where T is the age of the blast wave. The magnetic field on 
this scale can grow by at most of the order of one e-fold over 
this time. On the other hand, because the field is amplified 
along its original direction, scales below the maximum scale 
can grow exponentially, in contrast to mere stretching, which 
would continually alter the scale as the field lines lengthen. 

The formal objective of the present work is to show that the 
parameter a, as defined in standard dynamo theory, is non- 
zero in the presence of cosmic ray streaming instabilities, and 
to estimate its value. This will imply a maximum amount of 
growth on any given scale over any given time interval. 

The partial pressure P cr = n cr TmiC 2 /3 in cosmic rays in 
any logarithmic interval of energy above energy E m \ n is about 
P cr - pu 2 /ln(P max /P min ) < pu 2 J\Q (Ellison & Eichler 
1985). Here, u s is the streaming velocity, p = riirrii is the 
plasma density, n» is the interstellar number density of pro- 
tons, rrii is the mass of a plasma ion, and T is the CR's Lorentz 
factor. Thus, the high energy CR number density is 
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TniiC 2 TniiC 2 ]n(E max / E min ) ' 



The cosmic ray current J cr = n cr e u s due to CRs within a 
given logarithmic interval of CR energy can then be factored 
into dimensionless parameters as follows: 



4tt _ 3 u s P cr eB 
J ~ c kB~ 2 c B 2 /8ir kYm iC 2 ' 



(1) 



We will show in the present study that J' is a key parameter 
that determines large-scale magnetic field amplification. The 
last factor eB / (kTrriiC 2 ) in Equation (1) is the deflection of 
the CR over its passage through one scale length, of the 
magnetic field. The derivation of the Bell instability that is 
caused by the CR current in plasma, assumes that the deflec- 
tion of the CR is small. For the firehose instability, the subject 
of a separate investigation, it can be larger than unity, as the 
CRs would be bound to the field lines. However, a situation 
can occur in which the CRs maintain a steady anisotropy, A, 
even if they scatter, and the choice for the streaming veloc- 
ity is then u s — Ac, and the shock velocity is assumed to be 
approximately the same. 

The factor u s /c is typically of the order of 1CT 2 or less, 
both for the Galaxy as a whole, where the anisotropy is limited 
by direct measurements, and for the cosmic rays accelerated 
by forward shocks of blast waves from supernovae remnants 
(SNR), where the CR precursor moves essentially at the shock 
velocity, typically 10~ 3 to 10~ 2 times c. Very young SNR 
can have somewhat higher shock velocities, but the CR accel- 



erated in them may then be more prone to adiabatic losses. 
Note that because we assume the streaming velocity to be the 
velocity of the shock, we may allow for the possibility that 
CRs scatter in our estimate of the CR current. But we still 
have to limit the time over which a given parcel of fluid is 
exposed to the CR current to be less than the shock crossing 
time - the time it takes the shock to cross the length of the CR 
precursor at that energy. 

The factor 3P cr /{B 2 /Air) is of the order of unity for the 
Galaxy as a whole if all relativistic CRs are included. In this 
case, the anisotropy must be taken to be at most 1CT 3 and 
u s < 10 _3 c. In the case of supernova blast waves, the fac- 
tor iP" / (B 2 / Att) can be as high as <~ u 2 /v\, where va is 
the Alfven speed. Thus, the CR pressure can be a significant 
fraction of the ram pressure pu 2 ; see, e.g., Ellison & Eichler 
(1985). 

Altogether, we can choose a plausible value for J of 
~ 3v? 8 /cv\ ln(S max /S m i n ) which, for supernova remnants 
(SNR), can range from 10 4 for young SNR (u s ~ c/10 ~ 
10 3 ?;a) to order unity for old ones. For the Galaxy as a whole, 
we are probably limited to 3P cr / (B 2 / 4tt) < 1. For colli- 
sionless shock waves generated in the interstellar medium by 
GRBs, the value of J can be enormous - of order 10 14 , and 
the a effect may be particularly effective in that context. This 
will be discussed in greater detail below, where relativistic ef- 
fects are included more carefully. 

The physical time over which the mechanism can operate 
on a given patch of upstream fluid in the case of a SNR is 
one expansion time, or R s /u s , where R s is the radius of the 
forward shock. As u s does not appear explicitly in the simu- 
lation parameters, we use the above expression J <~ v? s /cv\ 
to substitute for u s in terms of J cr and va- For the Galaxy 
as a whole, we are limited by the Hubble time to be about 
10 4 _Rg/ w a, where Rq is the size of the Galaxy. 

2. GOVERNING EQUATIONS 

We consider magnetohydrodynamic (MHD) flows consist- 
ing of background plasma ions of number density rii, elec- 
trons of number density n e , and cosmic ray protons of num- 
ber density n cr . The equation of motion in MHD flows with 
cosmic rays imbedded in a background plasma, reads (Chen 
2010): 

P^r =-VP+^JxB + e(m-n e )E + F v , (2) 

where D/Dt = d/dt + U- V is the advective derivative, U is 
the fluid velocity, p and P are the fluid density and pressure, 
respectively, E and B are the electric and magnetic fields, and 
F v is the viscous force. The densities of the plasma current, 
J, and of cosmic ray protons, J cr , are the sources of magnetic 
field in Maxwell's equations: 



VxB = -(J + J cr ). 



Vx£ = - 



1 dB 



(3) 
(4) 



with Ohm's law: J = a(E + c~ 1 UxB), where a is the 
electrical conductivity of the gas, and we have neglected in 
Equation (3) the displacement current, because the conduc- 
tivity is high and the fluid motions slow compared with the 
speed of light. These equations yield the induction equation, 



— = VX (UxB - V VxB + cJ a ja) , 



(5) 
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where i] = c 2 /4na is the magnetic diffusivity. We assume 
quasi-neutrality for the whole system, i.e., rij + n cr = n e , and 
Equation (2) reads: 



p^- = -VP+ -i- (Vx-B) xB + F„- ij cr xB 



Dt 4tt 

+ ien cr (UxB), 
c 



(6) 



where we have used Equation (3) and assumed that | J/<r\ <C 
c~ 1 \UxB\, i.e., we consider plasma flows with large hydro- 
dynamic and magnetic Reynolds numbers. Hereafter we as- 
sume that the cosmic ray velocity is much larger than fluid 
velocity U, so that the last term in Equation (6) vanishes. The 
plasma density is determined by the continuity equation: 



dp 
dt 



+ V-pU = 0. 



(7) 



Following Bell (2004), let us consider the equilibrium: 

jeq + jcr = _gcq = B =CQm ^ [/cq = q an( j p cq = 

p» =const. Linearized Eqs. (3)-(7) for small perturbations 
yield the following dispersion relation: 



k 2 



■ u, 
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(8) 



where uja = k-VA is the frequency of the Alfven waves, 
uj s = kc s is the frequency of the sound waves, c s is the 

1/2 

sound speed, k is the wavenumber, va = £?*/ (47rp) ' is 
the Alfven speed, lu ci = c _1 J cr (Air/p) 1 / 2 , 7 = j B + V^ 2 , 
7b is the growth rate of an instability and we have considered 
the case where the equilibrium magnetic field B„ and J cr are 
directed along the z axis. 

In this system, the non-resonant Bell instability (Bell 2004; 
Lucek & Bell 2000) is excited by the cosmic ray current that 
causes growing MHD modes. The growth rate of this insta- 
bility for p > 1 that follows from Equation (8), is given by 



7b 



Ma w 



-1 1/2 



k 



uj 2 a 



T]k 2 



(9) 



where (3 is the ratio of gas pressure to magnetic pressure. For 
incompressible MHD modes b(k) — i(k-B^) u(k)/-/B (Bell 
2004), where u and b are perturbations of velocity and mag- 
netic field. 

When (3 <C 1, the growth rate of this instability that follows 
from Equation (8), is given by 



IB : 
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< l + V 2 
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- V k 2 



(10) 



The nonlinear stage of this instability has been investigated 
in a number of publications (see, e.g., Bell 2004, 2005; Pel- 
letier 2006; Reville 2008; Zirakashvili et al. 2008; Amato & 
Blasi 2009; Luo & Melrose 2009; Vladimirov et al. 2009; 
Zweibel & Everett 2010; Bykov et al. 201 la,b). These studies 
have demonstrated that this instability produces small-scale 



turbulence with cascading of turbulence energy into larger and 
smaller scales. 

In this paper we discuss the possibility of mean-field dy- 
namo action caused by the interaction of mean electric cur- 
rent of cosmic ray particles with small-scale background ho- 
mogeneous turbulence produced by the Bell instability. This 
paper can be considered as an extension of the recent study 
by Bykov et al. (2011a) who demonstrated, using a multi- 
scale quasi-linear mean-field approach, that small-scale Bell- 
type turbulence can result in the growth of long-wavelength 
obliquely propagating modes (Bykov et al. 201 la). 

3. LARGE-SCALE INSTABILITY 

In this section we discuss mean-field dynamo action in 
small-scale turbulence produced by the Bell instability in a 
plasma with a given mean electric current of cosmic ray ions. 
The importance of this effect is determined by the ratio of the 
cosmic ray current to the ambient field; see Equation (1). 

3.1. Mean field dynamo equations 

We use a mean field approach in which magnetic and ve- 
locity fields are divided into mean and fluctuating parts, U = 
U + u and B = B + b, where u and b are fluctuations 
of velocity and magnetic field, B and U are the mean mag- 
netic and velocity fields, and the fluctuating fields have zero 
mean values. We consider the case when magnetic and fluid 
Reynolds numbers are large. This implies that the nonlinear 
terms in the induction and Navier-Stokes equations are much 
larger than the dissipative and viscous terms. In this case the 
quasi-linear approach for determining the turbulent transport 
coefficients (e.g., the a effect and turbulent magnetic diffusiv- 
ity) does not work. We use instead the spectral r relaxation 
approximation (Orszag 1970; Pouquet et al. 1976; Kleeorin 
et al. 1990; Rogachevskii & Kleeorin 2004) that is valid for 
large magnetic and fluid Reynolds numbers. A justification 
for the t approximation in different situations has been per- 
formed in numerical simulations and analytical studies (see, 
e.g., Brandenburg & Subramanian 2005; Rogachevskii et al. 
201 1). For more details see Appendix A. 3. 

Averaging Eqs. (5) and (6) over an ensemble of turbulent 
eddies yields the following mean-field equations: 



dB 
~dt_ 
DU 

157 



= Vx (UxB + uxb-rjVxB) , 



T> = -VP + — (VxB) xB- -J"xB 

47T v ' c 



(UxB) - VjTZuJ + Fv, 



(11) 



(12) 



where S(B) = u x b is the mean electromotive force, «7 cr 
is the mean density of the electric current of cosmic ray par- 
ticles. For large hydrodynamic and magnetic Reynolds num- 
bers we can neglect kinematic viscosity, v, and magnetic dif- 
fusivity, 77, in comparison with the turbulent viscosity and tur- 
bulent magnetic diffusion. 

3.2. Contributions to the a effect 

We will show in this study that, formally, there are two 
contributions to the a effect caused by: (i) existing ki- 
netic helicity produced by the Bell-instability (referred to as 

tt(°) • (V x below; (ii) correlations in the forcing by 

the mean cosmic ray current in the presence of small-scale 
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anisotropic magnetic fluctuations which create further pertur- 
bation k/ 1 ) in the velocity field u^°\ Effect (ii) causes oppo- 
site sides of a magnetic loop to be forced in opposite direc- 
tions, thereby twisting the loop out of its original plane. The 
distinction between the two terms may be more formal than 
physical. That is, in an unstable circularly polarized Alfven 
mode, the additional helicity added by the cosmic ray forcing 
term is merely a continuation of the process that formed the 
already existing helicity. 

The a effect, however, is distinct from the linear growth 
of the circularly polarized Alfven wave. It can be thought 
of as one circularly polarized Alfven wave riding on another, 
somewhat longer wavelength Alfven wave (which need not 
be circularly polarized). The longer wave makes a perpen- 
dicular component to the original field, while the stretching 
of the perpendicular component (into a loop, say) together 
with the twisting of the loop restores newly created flux back 
into the original direction. It is the nonlinear coupling of two 
waves, each of which is of the sort discussed by Bykov et al. 
(201 la). This coupling can scatter energy in the two modes of 
wavenumbers fci, fc 2 into modes of much longer wavelength, 
and thereby amplify the field at large scale in its original di- 
rection. 

Let us consider the case where the equilibrium uniform 
mean magnetic field B» and the mean density of the elec- 
tric current of accelerated particles J cr are directed along the 
z axis. We take into account effects which are linear in per- 
turbations of the mean magnetic field: B — B — JB*, i.e., we 
consider kinematic mean-field dynamo. 

The first contribution to the a effect is caused by the helical 
part of the turbulence. A non-zero kinetic helicity is caused by 
the Bell instability that results in the production of small-scale 
helical turbulence. This contribution to the mean electromo- 
tive force is given by £> — a\) Bj + ... (see Appendix A.6), 



where af) = af Sij is an isotropic a effect, 



a5 r = -CiT u(°)-(Vx«(°)) ) 



(13) 



with coefficient C\ w 1/3, and dots referring to higher-order 
terms that will be considered later in Section 4.2. In the rest of 
this section the dots will not be noted explicitly. Here, are 
the velocity fluctuations of the background turbulence. Equa- 
tion (13) is a well-known result for the a effect caused by 
kinetic helicity of the turbulence (see Krause & Radler 1980; 
Moffatt 1978; Parker 1979; Zeldovich et al. 1983). 

The second contribution to the a effect is caused by incom- 
pressible and non-helical parts of the anisotropic turbulence 
interacting with the mean cosmic ray current. In particular, 
this contribution to the mean electromotive force is given by 



Bj, where afp = af (Sij + e^j) (see Ap- 



pendix A.6), 



a c J = C 2 



4tt J c 



-B* 



1/2 



V A sgn ( J" • B» 



(14) 



V A = B^/iAnp) 1 / 2 is the mean Alfven speed based on the 
equilibrium mean magnetic field B*, e is the unit vector di- 
rected along £?«, C 2 = 4e (q - l)/3(2g - 3), < e < 1 
is the anisotropy parameter of small-scale turbulence and q is 
the exponent of the energy spectrum of the turbulence. Note 
that C 2 w 8/3 for e w 1 and q = 5/3. The total a effect is 



(15) 



to the a 



the sum of the two contributions: 

a cr = a (i) + J U) 

v »j »j 

The mechanism of the second contribution, 

"j 

effect can be understood by the following reasoning. Tangling 
of the mean magnetic field B = B — B* by velocity fluctua- 
tions of the background anisotropic turbulence produces 
magnetic fluctuations: 

^«(B.V)u<°>. (16) 

The generated magnetic fluctuations &W interacting with the 
cosmic ray current J cr produce additional velocity fluctua- 
tions: 



— - oc = J cr xb (1) . 

at cp 



(17) 



These velocity fluctuations contribute to the mean electro- 
motive force £^ IV> = itW x &(°). Here, are the mag- 
netic fluctuations resulting directly from the Bell instability, 
just like the velocity perturbations u^°\ In particular, the 
stretching of the original magnetic field determined by Equa- 
tion (16), and rotation of the stretched magnetic loop deter- 
mined by Equation (17) create an electric field £^ IV> along the 
original magnetic field B. The Lorentz force in Equation (17) 
plays the role of rotation of a magnetic loop. Using Eqs. (16) 
and (17), we can estimate the mean electromotive force £^ u ' 
using dimensional reasoning. Indeed, the velocity fluctuations 
can be estimated as 

2 

«W _^_J^ xb (i) K _I_ j^x(B-V)u {0) , (18) 

cp cp 

where r is the characteristic time of the turbulence. Therefore, 
the mean electromotive force £^ is estimated as 



^ (II) oc 



■ 3"i 



) Bj. (19) 



The mean electromotive force £^ is proportional to the 
mean magnetic field B, i.e., £ 



(ii) 



af^ is the symmetric part of . Therefore, 



a,j Bj and the tensor 



_75? 6( )v - u (0) - J 



(°)y7.„.(°; 

3 



(20) 



It follows from Equation (20) that for isotropic background 
turbulence, b^u-p oc 6{j, and afp vanishes. Note also that 

in many kinds of background turbulence, the tensor b^u^ 

as well as the background mean electromotive force £^ = 
M (o) x fo(°) vanish. On the other hand, for the turbulence pro- 
duced by the Bell instability, the tensor b^u^ does not van- 
ish. 

Now let us compare this mechanism for the a effect with 
the standard mechanism based on nonzero kinetic helicity. 

The mean electromotive force £^ = «(°) x bW, where the 
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magnetic fluctuations are generated by the tangling of 
the mean magnetic field B by the velocity fluctuations of the 
background turbulence «/ '. Therefore, the mean electromo- 
tive force is given by 



f (I) ar«(°) x (B-V)u(°), 



(21) 



and the o$ tensor is 



- /// in 

u^ViU^ ), (22) 



where the kinetic helicity is produced by the Bell instabil- 
ity and Xij = EmmllmV ' 3 Un^ + £jnm UmV iUn^ is the 

symmetric tensor of the kinetic helicity for anisotropic tur- 
bulence. Using the model for helical background turbulence 
Equation (22) can be reduced to Equation (13). Therefore the 

afj tensor is directly related to the kinetic helicity. Now let 
us estimate the kinetic helicity for the Bell background tur- 
bulence. The contribution of the cosmic ray current to the 
velocity of the background turbulence is 



duW 1 -T- u(0) 

— — - oc — -J cr x6 (0) . 
at cp 



(23) 



Therefore, the corresponding contribution of the cosmic ray 
current to the vorticity of the background turbulence is 



^Vx/'oci (J^-V)&<°>. 

dt cp v ' 



(24) 



Multiplying Equation (23) by the vorticity Vxti' ' and Equa- 
tion (24) by the velocity u^°\ and averaging over the ensem- 
ble we obtain: 



d- 

0? 



JV- 

(V x «(0)) oc --^ x (V x 



cp 

-Un Vj 



(25) 



Dimensional reasoning yields the following estimate for the 
the kinetic helicity for the Bell background turbulence: 



«(°)- (V x u(°)) oc 



J? 



C 2 p 2 



e 6 (0) V-6 



(0) 



+- 



cp 



u (0) V-6 (0) 



(26) 



where we h ave taken into account that V x oc 
(r/c~p) (J CI •V)fo(°- 1 . The first term in the right hand side of 
Equation (26) is related to the current and magnetic helici- 
ties. For the Bell instability &(°)(fc) w i(k-B«,) «(°) {k)/ lB . 
Therefore, Eqs. (20) and (26) show that only part of the con- 
tribution to afp may be related to the kinetic helicity. This 
contribution is caused by the last term in the right hand side of 
Equation (26). Note that Equation (26) for the kinetic helicity 
is different from Equation (C6) of Bykov et al. (201 la), where 
no substitution for b^\ as in Equation (18), is performed. 

3.3. Large-scale instability for B(t, x, y) 

Let us start to analyze the large-scale instability using a 
case of incompressible flow when perturbations of velocity 



and magnetic field are independent of z. In this case pertur- 
bations of the mean magnetic field B(t, x 7 y) are determined 
by the following equation: 



a r> 

— = (b*-v)u + vx ni^b) 

ot v ' 



(27) 



where (uxb). = g. _ ^(VxJB);, i] t = C^u ^o is 
the turbulent magnetic diffusivity, C v ~ 1/3 is a constant, 
lo is the maximum (integral) scale of turbulent motions, u 
is the characteristic turbulent velocity in the integral scale of 
turbulence, B = B*+ B and U(t, x, y) are the perturbations 
of the mean velocity field. Here, for simplicity, we neglect 
small anisotropic contributions to the turbulent magnetic dif- 
fusion. Since the equilibrium uniform mean magnetic field 
B* is directed along z axis, the first term, (£?»-V)t/, on the 
right hand side of Equation (27) vanishes for perturbations 
which are independent of z. In this case the mean-field in- 
duction equation (27) is decoupled from the mean-field mo- 
mentum equation. The perturbations of the mean magnetic 
field B can be written in the form of the axisymmetric field: 
B = B z (t, x,y)e + Vx[A{t, x, y) e}. 

Let us start the analysis with the simpler case in which 
the tensor af^ = a cr Sij is isotropic. Then the functions 

B z (t,x,y) and A(t,x,y) are determined by the following 
equations which follow from Equation (27): 



dB z 
dt 
OA 
~dt 



= -a cr Ai + 77 f A n B Zl 
= a cr B z +T lt A H i, 



(28) 
(29) 



where Ah = V 2 + V 2 . We seek a solution of the mean-field 
dynamo equations (28) and (29) of the form cx exp[7i nst t + 
i(K x x + K y y)\. The growth rate 7 inst of the mean-field dy- 
namo instability in homogeneous turbulent plasma with cos- 
mic rays is then given by 



7inst = \a"K\- m K 2 



(30) 



where K 2 = K 2 + K 2 and a cr = a". This large-scale dy- 
namo instability is called a 2 dynamo, because this dynamo 
is caused by the interaction between two modes, toroidal, 
B^ — B z (t,x,y) e, and poloidal mean magnetic fields, 
where poloidal field, B^ = Vx[i(t, x, y) e] is determined 
only by the potential A(t, x, y). The toroidal field is gener- 
ated from the poloidal field by the a effect due to the first 
term, — a cr AA in Equation (28), while the poloidal field is 
generated from the toroidal field by the a effect [due to the 
first term, a cr AB z in Equation (29)]. This implies that the a 
effect acts twice in the positive feedback loop, and these inter- 
actions between the magnetic field components cause the a 2 
mean-field dynamo. 
Now we consider the case in which the tensor afi is 

anisotropic. In this case the functions B z (t,x,y) and 
A(t, x, y) are given by: 



dB z 
dt 



= -a y I y V 2 x A-a x r x V 2 y A + r 1t A H B z , (31) 

(32) 



dA 

— =af z B z +r] t A H A 
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where a" x = a" y ^ af z . The growth rate 7i nst of the mean- 
field dynamo instability in this case is 



7inst = \K\ ^"yy a ?z ~ Vt K , 



(33) 



and the ratio of magnetic energies along and perpendicular to 
the direction of the cosmic ray current is given by 



Bl 



B 2 



B 2 

v 



(34) 



where 



*yy — "^1 ~ r — u i "i" ■ 

Comparing with Bykov et al. (201 la), who studied a long- 
wave instability for incompressible flows of a plasma with 
cosmic rays, we note that their growth rate vanishes when 
perturbations of mean magnetic and velocity fields are inde- 
pendent of z; see their Equation (23). In this sense, the mean- 
field dynamo mechanism studied in our paper for K z = is 
a complementary effect to a mechanism related to the long- 
wave instability discussed by Bykov et al. (201 la). 

3.4. Large-scale instability for B(t, z) 

Now let us consider the case when perturbations of mean 
magnetic and velocity fields depend only on z. The pertur- 
bations B z = (because V • B = 0), and the large-scale 
dynamo instability can generate magnetic field only in the di- 
rection perpendicular to the cosmic ray current. In this case 
we seek a solution of the linearized mean-field dynamo equa- 
tions (1 1) and (12) of the form oc exp(7; nst t + iK z z). The 
growth rate 7; nst of the large-scale dynamo instability is given 
by 



2a% 



Tinst — 



\n A n cr \-n 2 A + ^( a "K z f 



1/2 



i 

+ 2 



a cr K z \- Vt K 2 z 



(35) 



where fL 



k z v a , n c 



cr 1 J cr (4.7T/P) 1 / 2 and a c 



a" = a C y y — af + a™. The growth rate 7i nst of the dynamo 
instability can be interpreted as the interaction of the a 2 large- 
scale dynamo instability [determined by the terms oc a CT K z 
in Equation (35)] and the Bell instability [described by the 
terms oc \£Ia ^ cr | — Q\ m Equation (35)] in a homogeneous 
turbulent plasma with a cosmic rays current. When a cr = 0, 
the growth rate of the large-scale instability (35) coincides 
with that derived by Bykov et al. (2011a); see Eq. (23) of 
their paper. 

4. DIRECT NUMERICAL SIMULATIONS 

4.1. DNS model 

We consider a cubic computational domain of size L 3 . The 
smallest wavenumber is fci = 2-k/L. We adopt an isother- 
mal equation of state with constant sound speed c s , so the gas 
pressure is p = pc 2 . We solve the equations of compressible 
MHD in the form 



DC/ _ 1 



(VxB)xB - c 2 Vp + V • (2upS) 



--J cr x B, 

c 

OA 

— = UxB + 7 1 V 2 A, 



(36) 
(37) 
(38) 



where v and rj are kinematic viscosity and magnetic diffusiv- 
ity, respectively, B = Bo + V x A is the magnetic field con- 
sisting of a uniform mean background field, B = (0, 0, B ), 
and a nonuniform part that is represented in terms of the mag- 
netic vector potential A, andSy = ^(Uij + Uj^i)— ^Sij'V -U 
is the traceless rate of strain tensor, where commas denote par- 
tial differentiation. 

In all cases we adopt triply periodic boundary conditions. 
The simulations are performed with the Pencil Code 1 , 
which uses sixth-order explicit finite differences in space and 
a third-order accurate time stepping method (Brandenburg & 
Dobler 2002). Simulations have been done with various reso- 
lutions, but here we focus on two runs with 512 3 mesh points. 
Some of the results are also compared with corresponding 
ones at 256 3 mesh points. 

4.2. Test-field method 

We apply the quasi-kinematic test-field method (see, e.g., 
Schrinner et al. 2005, 2007; Brandenburg et al. 2008) to de- 
termine all relevant components of the tensor ay and turbu- 
lent magnetic diffusion. This method allows for the presence 
of strong magnetic field as long as the magnetic fluctuations 
are entirely a consequence of the imposed field (Rheinhardt 
& Brandenburg 2010). The essence of this method is that a 
set of prescribed test fields B^ p ' q ^ and the flow from the DNS 
are used to evolve separate realizations of small-scale fields 
b<P'i\ Neither the test fields nor the small-scale fields 

lf(p,q) ac t back on the flow. These small-scale fields are then 
used to compute the electromotive force £(p^ corresponding 
to the test field B^ p ' q \ The number and form of the test fields 
used depends on the problem at hand. 

The choice of test fields depends on the averaging that is 
performed. Relevant for the present study are averages that 
depend on x or y, or both. To gather sufficient statistics, we 
adopt planar yz averages that depend on x, so we use test 
fields £?( lc ) = {0,B coskx,0),W ls ^ = (0, B sin kx, 0), 
I?( 2c ) = (0, 0, B cos kx), andl?( 2s ) = (0, 0, B ,sin kx), 
in which case the series expansion of the electromotive force 
contains two terms 



ciij Bj 



t]ij{V x B)j. 



(39) 



The symmetric part of the tensor is of particular interest 
and is commonly referred to as the a tensor, 



a. 



— , {o-ij + a j 



(40) 



Errors are estimated by dividing the time series into three 
equally long parts and computing time averages for each of 
them. The largest departure from the time average computed 
over the entire time series represents an estimate of the error. 

4.3. DNS results 

All runs are isothermal with c s = 10. The box has a size 
of 167r, i.e., ki — 1/8. In our units, An/c = 1, po = 1 and 
the magnetic Prandtl number v/tj = 1. The relevant non- 
dimensional parameters are Lu = vao/v^i ( tne Lundquist 
number), where «ao — is the non-dimensional Alfven 
speed, and J = 4ir J cr / ' cBok\ (the non-dimensional cosmic 
ray current density). 

Visualizations of the magnetic field B x /B on the periph- 
ery of the computational domain are shown in Figure 1 for 

1 http : //pencil-code . googlecode . com 
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FIG. 1 . — Visualization of B x /Bo on the periphery of the computational domain using 512 3 mesh points for J cr = 0.1, Bq = 0.01, fcj = 1/8 (so that 
J = 80), and v = r) = 10~ 3 (so that the Lundquist number Lu = 80). 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 

t v k0 k l t v k0 k l 



FIG. 2. — Time evolution of E^ki/v^ for modes with different 
wavenumbers for the run with J = 80. The short straight lines show the 
growth rate of the Bell instability, as given by Equation (9) for modes with 
three selected values of k, as well as the value of 7i nst as given by Equa- 
tion (33). 

a run with 512 3 mesh points using the following parameters: 
the non-dimensional cosmic ray current density is J = 80 
and the Lundquist number Lu = 80 (with J cl = 0.1, Bq = 
0.01 and fc x = 1/8). In the beginning of the instability, the 
length scale is rather small, but it increases continuously as 
time goes on. Note in particular the much larger horizontal 
length scales in the xy plane that may be associated with the 
dynamo instability. After tvAoki « 0.5 the instability reaches 
yet another stage during which the magnetic field pattern be- 
comes turbulent. Again, as time goes on the typical eddy scale 



FIG. 3. — Time evolution of the ratio of the spectral vertical (along the 

imposed field Bo) and horizontal magnetic energies 2E^/E^ for the ran 
with J = 80. 

increases. 

All runs with a constant cosmic ray current show that 
there is a growth until the velocities become eventually su- 
personic. This is probably the reason the simulation termi- 
nates. It is conceivable that this could be avoided by includ- 
ing the backreaction of the amplified field on the cosmic ray 
current, which would limit the Bell instability (Riquelme & 
Spitkovsky 2009). In Figs. 2 and 3 we show the time evo- 
lution of the normalized spectral vertical magnetic energies 

e\ k\ I v\ and the ratio of the spectral vertical (along the im- 
posed field B ) and horizontal magnetic energies 2e\j/E^ 
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for modes with different wavenumbers for the same run as in 
Figure 1 . As follows from these figures the dynamics of the 
instability has the following stages: 

(i) In the early stage there is the development of small- 
scale instability that results in the production of small- 
scale turbulence. It is seen in Figure 2 that the mode- 
averaged growth rate of this instability in this stage is 
slightly smaller than that of Equation (9) that describes 
the Bell instability for the fastest growing mode of a 
given \k\, as should be expected. 

(ii) In the second stage, there is formation of large-scale 
magnetic structures (during the interval 0.2-0.4 Alfven 
times; see Figure 1). Note from Figure 2 that the growth 
rate of the large-scale mode k = ki is about twice the 
growth rate of the largest mode. The interpretation of 
this will be given below. 

(iii) In the final stage, there is a development of larger- 
scale turbulence; the perturbed field actually exceeds 
the original field by a considerable factor, and a sig- 
nificant fraction of the energy is now present in modes 
with k± ~ A: 1 1 . The growth rate in this final stage agrees 
with that predicted by Equation (33), where a y r y , a c zz , 
and r\ t have been obtained with the test-field method, as 
described below. 

Time evolution of power spectra of magnetic energy of the B z 

component, E^, and those of all components, Em, are shown 
in Figure 4, which demonstrates an inverse energy cascade- 
like behavior. The evolution during the second stage (during 
the interval 0.2-0.4 Alfven times; see Figure 1), shows that 
large scale modes (k ~ kx) are growing at about twice the 
growth rate of the fastest growing mode k/ki ~ 40). We in- 
terpret this as perturbation field growth db/dt at large scale 
due to the coupling of pairs of higher k modes (fci and k^, 
say), and therefore proportional to M^ ms i3o oc exp(7i + 72), 
where 71 and 72 are the respective linear growth rates of 
the high k modes. During this stage, the perturbed field is 
still small compared to the original field, so the exponential 
growth of the total field as described by equations (13) and 
(14) has not yet begun. 

At the end of the simulation, both kinetic and magnetic en- 
ergy spectra develop a fc~ 5 / 3 energy spectrum. This is shown 
in Figure 5, where we plot spectra compensated by e~ 2 / 3 fc 5 / 3 , 
where e is the total (kinetic and magnetic) energy dissipa- 
tion rate. The wavenumber is normalized by the dissipation 
wavenumber kd = [e/ {y + r;) 3 ] 1 ' 4 . 

The test-field results for the normalized components of the 
a tensor, a zz vao /u 2 ms and a yy vao /u 2 ms are shown in Fig- 
ure 6 for the same run. The value of a™ s VAo/u 2 ms « 0.6 is 
of the same order of magnitude as that determined by Equa- 
tion (14), i.e., a\ l z coly vao /n, 2 ms ~ 0.5. Note that at late times, 
after about 0.45 Alfven times, we find that both a™ /u Ims 

and a™ /« rms are approximately constant in time. This 
measured value of the a effect is much larger than that based 
on kinetic helicity determined by Equation (13), unless tq is 
considerably larger than l/ku Tms see Figure 7. We interpret 
this as evidence for an additional contribution, Equation (14). 
The growth rate of the large-scale dynamo instability is of the 
same order as that theoretically predicted from the large-scale 
dynamo instability. Indeed, Equation (33) for the growth rate 
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FIG. 4. — Time evolution of E M and Em for the run using 512 3 mesh 
points (J cr = 0.1, Bo = 10~ 2 , so that J = 80). The normalized time 
interval between different spectra is v\okiAt = 3/80 ~ 0.04. The solid 
lines refer to the initial spectra proportional to fc 4 for small values of k and 
the red and blue lines represent the last instant of Em and E^ , respectively. 
The straight lines show the k A and fc~ 5 / 3 power laws. 
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FIG. 5. — Compensated spectra of Ek (blue), Em (red), and E' M (black), 
at the end of the simulation. Here, e is the total (kinetic and magnetic) energy 
dissipation rate. The dashed horizontal line goes through 2.7. 

yields 7 thoor y w 0.03, corresponding to 7 thcory /wac^i ~ 24 
in non-dimensional units. This value is in agreement with the 
growth rate expected from an a 2 dynamo with coefficients 
obtained with the test-field method; see the inset of Figure 7, 
where 7i n st/vAofci ~ 20. This agrees with the growth seen 
in the DNS at k = k\ during the time interval 0.45-0.65; 
see Figure 2. Finally, Equation (34) for the ratio of magnetic 
energies along and perpendicular to the direction of the cos- 
mic ray current yields [B 2 z /(Bl + B 2 )] thc ° ly 0.08, while 
this ratio according to DNS (see Figure 3) is estimated as 

[B 2 z /(Bl + B 2 )] UNS « 0.06. The above comparison im- 
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FIG. 6. — Normalized ttzz^Ao/^Ls an< ^ a yy »Ao/ u nns f° r me mn 
with J = 80. The inset shows that after about 0.45 Alfven times, both 
a H > y ,S / nrms an ^ "If 8 /"'" are approximately constant in time. 
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FIG. 7. — Evolution of a = (a yy a zz ) 1 ^ 2 , normalized by t ci (lj-u) where 
r cr = l/u) CT for the run with J = 80. The inset gives the instantaneous 
value of 7i nst as derived from Equation (33). 
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FIG. 8. — Same as the bottom panel of Figure 4 (J = 80), but showing 
only the time interval 0.35 < tv^oki < 0.45 (dashed lines), i.e., the last 
0.1 Alfven times just near the end of the linear growth phase. The black solid 
lines refer to tv^ki = 0.35, while red lines and blue lines refer to kf(t) 
and fcjr (t) at tv^ki = 0.65. The slope fc 3 is shown for comparison. 

plies a good agreement between theoretical predictions and 
the results of the numerical simulations. 
In Figure 8 we show the time evolution of power spectra of 

magnetic energy of the B z component, E M , and those of all 
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FIG. 9.— Evolution of kf (solid, red), kf^ (dashed), and kf (dotted, 
blue) for the same run as in Figure 1 ( J = 80). 

components, Em, like it was presented in Figure 4, but during 
only the last 1/10 of an Alfven time. We see that the ampli- 
fication of magnetic field with respect to initial perturbations 
during the last 1/10 of an Alfven time (which is of the order 
of the shock crossing time) at k = k\ is 3.5 orders of magni- 
tude (i.e., up to 0.6 x 10~ 3 of the field). However this field 
is less than the equilibrium imposed field Bq = 10~ 2 . On the 
other hand, the visualizations of the magnetic field B x /B on 
the periphery of the computational domain shown in Figure 1 , 
demonstrates that the ratio B x /B < 60 at tt>Ao^i — 0.65. 

To quantify the inverse cascade-like behavior, let us now 
look at the evolution of the wavenumber of the magnetic 
energy-carrying eddies, kf, defined via 

[/c f M (i)] _1 = J k- 1 E M (k,t)dk/ J E M (k,t)dk, (41) 

and likewise for kinetic energy as well as magnetic energy in 

the z component, kf and kf* , respectively (see Figure 9). It 
turns out that all three wavenumbers reach a value somewhat 
above 20 k% by the end of the small-scale dynamo instability, 
and then drop rapidly in the mean-field dynamo stage. Note, 
however, that the decrease of kf is somewhat faster than that 
of/c f M . 

The test-field method yields not only ay, but also the com- 
ponents of the turbulent magnetic diffusivity tensor rjij. It 
turns out that r\ vy w r] zz . The ratio %/r], where r\ t — 
(rjyy + i] zz )/2, greatly exceeds unity in the late stages; see 
Figure 10. Mean-field dynamo efficiency depends on the dy- 
namo number D = (a/rjTki) 2 , where a = {ayyazz) 1 / 2 and 
Vt = Vt + V- The dynamo number exceeds the critical value 
for dynamo action of unity (D > 1) after about 0.3 Alfven 
times; see the inset of Figure 10. 

4.4. Comparison with the run J — 800 

In this subsection we discuss the results of DNS with the 
higher cosmic ray current, i.e., when the normalized cosmic 
ray current is increased by one order of magnitude J = 800 
(i.e., J cr = 1, B = 0.01 and fc x = 1/8), and the other 
parameters are the same as in the previous subsection (the 
Lundquist number Lu = 80 and the resolution is 512 3 mesh 
points). In this case very small scales (higher k modes) are 
not well resolved. On the other hand, the main contribution 
to the mean-field dynamo (which is the main subject of our 
study) is of the maximum scale of turbulent motions. The 
time evolution of spectra of magnetic, Em, and kinetic, Ek , 
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FIG. 10. — Evolution of r]t/r] (solid line), rjyy/r) (dotted line), and rj zz /r] 
(dashed line) for the run J = 80. The inset shows a/rjtki. 
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FIG. 13.— Evolution of B rms / p 1//2 E/ r ms for runs with — 80 (dotted 
and dashed lines for 256 3 and 512 3 meshpoints, respectively) and J = 800 
(solid line, 512 3 meshpoints). Upper horizontal axis corresponds to case 
J = 80, while lower horizontal axis corresponds to case J = 800. 
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FIG. 1 1 .— Time evolution of E M and E K for the run J = 800 (J" = 1, 
B = 10~ 2 ) showing only the time interval 0.0125 < tv A0 ki < 0.075. 
The time interval between the lines At^Ao^l = 0.0125. The red and blue 
lines represent the last instant of Em and Ek > respectively. The slopes fc 4 
and A; 3 are shown for comparison. 
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FIG. 14. — Instantaneous growth rate 7 DNS = d In C/ rms /dt (upper panel) 
and ln[/rms/ffci (lower panel) for runs with S = 80 (dotted and dashed 
lines for 256 3 and 512 3 meshpoints, respectively) and J = 800 (solid line, 
512 3 meshpoints). Upper horizontal axis corresponds to case J = 80, while 
lower horizontal axis corresponds to case J = 800. 



FIG. 12. — Normalized B rms /BQ for runs with J = 80 (using 256 3 and 
512 3 meshpoints; dotted and dashed lines, respectively) and J = 800 (512 3 
meshpoints, solid line). Upper horizontal axis corresponds to case J = 80, 
while lower horizontal axis corresponds to case J = 800. 



energies for J — 800 is shown in Figure 11, that demonstrate 
an inverse energy cascade-like behavior, but to a lesser extent 
in comparison with the lower CR current ( J = 80). 



For comparison with the case J — 80, in Figs. 12-14 we 
show the time evolution of the total magnetic field B rms /Bo, 
^rms/p^^rms (which does not include the imposed field 
Bo) and the growth rate of the total velocity field, ^< DNS = 
d In U lms / dt, for runs with J — 80 (dotted and dashed lines 
for 256 3 and 512 3 meshpoints, respectively) and J = 800 
(solid line). Since the time evolution for different values of 
the cosmic ray current occurs on different time scales, we use 
the upper horizontal axis for the case J = 80, while the lower 
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FIG. 15. — Normalized a zz /U r ms for runs with — 80 (dotted and 
dashed lines for 256 3 and 512 3 meshpoints, respectively) and J = 800 
(solid line, 512 3 meshpoints). Upper horizontal axis corresponds to case 
J = 80, while lower horizontal axis corresponds to case J = 800. 

horizontal axis corresponds to case J = 800. 

Inspection of Figure 12 shows that, at the final stage of evo- 
lution, the generated magnetic field is by one order of magni- 
tude larger than the imposed field B . The J = 80 results at 
lower resolution are similar to those at higher, but have run for 
a slightly shorter time before resolution problems occurred. 
The growth rate of the velocity and magnetic field is 5 times 
larger for the case of J — 800 (see the upper panel of Fig- 
ure 14). On the other hand, the evolution of the kinetic energy 
- apart from this rescaling of time - is not strongly dependent 
on the cosmic ray current (see the lower panel of Figure 14), 
and neither is a zz /U rms (see the test-field results of the mea- 
sured a effect in Figure 15 shown in log scale). 

As noted above, the quasi-kinematic test-field method is 
valid as long as the magnetic fluctuations are entirely a con- 
sequence of the imposed field. In one particular case we have 
verified this by comparing with results from a fully nonlinear 
test-field method where velocity fluctuations resulting from 
the interaction with magnetic field fluctuations are also in- 
cluded. This method has currently been tested and imple- 
mented in the PENCIL CODE only for a modified set of equa- 
tions in which the pressure gradient and the U ■ VU term 
are omitted, but the Lorentz force is fully retained (Rhein- 
hardt & Brandenburg 2010). We have applied this method, 
with the modified set of equations, to a case similar to that 
displayed in Figure 6, but at lower resolution (64 3 ). In that 
case, a yy VA/U 2 ms is nearly constant after tvAki =0.1 and 
comparable to the corresponding value shown in Figure 6 at 
tvp^ki = 0.4, while a zz VA/U 2 ms agrees with that of Fig- 
ure 6 in the full time interval. More importantly, however, the 
quasi-kinematic and fully nonlinear test-field methods agree 
with each other within machine precision, confirming thus 
the applicability of the quasi-kinematic method to the present 
case. 

4.5. Interpretation of the Results 

Our results seem consistent with the following interpreta- 
tion. There appear to be three distinct stages. In the first stage, 
the c~ 1 J CI x b force (which we refer to as the Lorentz force 
due to the "counter-CR current" - i.e. the current in the ther- 
mal plasma that cancels the cosmic ray current) amplifies the 
motion perpendicular to the original, unperturbed magnetic 
field in one circular polarization of Alfven modes, thereby 



stretching the field such that it develops a component that is 
perpendicular to the original direction. This is the Bell insta- 
bility in our simulation, but anything that creates a perpendic- 
ular component might work just as well for this phase. The 
Bell instability grows fastest on small scales and in a direction 
whose k vector is parallel to the z axis. Tentatively, we may 
attribute the increasing preference for "perpendicular" energy 
that we see in the simulation, i.e. energy in motion in the "per- 
pendicular" (x and y) directions, to the faster growth rate of 
on-axis (k = k z ) Alfven waves. However, there is clearly 
significant parallel energy, even in the linear growth regime, 
presumably due to off-axis waves. The reason that the ratio 
of "parallel" to "perpendicular" energy decreases in the first 
stage is presumably because the on-axis waves, which have 
only motion perpendicular to the axis, are the fastest growing 
modes. 

During the first stage, the ratio of a zz to u rms grows expo- 
nentially, because for any given mode, aa«'Vx«~ ku 2 
is proportional to it 2 , and u is growing exponentially (Fig- 
ure 14), consistent with the linear Bell instability. The cor- 
relation length of both the magnetic field and the turbulence, 
apart from a drop in the very early stages, remains more or 
less constant, and corresponds to the scale of the fastest grow- 
ing modes. 

The second stage is much like the first, except that the 
growth rate of the larger scale, low k, modes begins to in- 
crease above its linear value; see Figure 2. This is presum- 
ably due to the fact that their growth is dominated by nonlin- 
ear coupling of higher k modes, which have developed much 
larger amplitudes than the low k ones, even though they are 
still in the linear regime. 

The third stage, for J = 80 (J = 800), begins at about 
tvAoki ~ 0.45 (tvAoki ~ 0.06), as shown in the upper panel 
of Figure 14. Several things clearly happen at the onset of the 
second stage: (i) The growth of U lms suddenly slows (Fig- 
ure 14), (ii) the growth of the magnetic field also shows a 
change (Figure 12), though it continues to rise, (iii) the cor- 
relation lengths of both the magnetic field and the velocity 
begin to increase noticeably (Figure 9), (iv) the ratio of a to 
w rms stops rising exponentially and either flattens out or grows 
much more slowly (Figs. 6 and 15), (v) the turbulent magnetic 
diffusivity begins to rise significantly and dominates the mi- 
croscopic value (Figure 10), (vi) the ratio of energy in parallel 
magnetic field to that in perpendicular field rises sharply (Fig- 
ure 3), and (vii) a Kolmogorov-type spectrum is reached from 
below, while the level of turbulence grows more slowly. 

All these changes can be understood in terms of nonlinear 
effects. Once the turbulent velocity exceeds the Alfven veloc- 
ity, the nonlinear convective term in the Navier-Stokes equa- 
tion becomes as or more important than the Lorentz force due 
to the counter-CR current, so that the stirring of the fluid by 
the latter, as expressed in unstable Alfven modes, is in equi- 
librium with eddy viscosity. By the same token, the amplitude 
of the magnetic field is large, so that the a 2 dynamo is acti- 
vated. Taken in isolation, an exponentially growing a effect 
would lead to super-exponential growth, but this is not what 
is seen. The sudden rise of parallel to perpendicular magnetic 
energy could perhaps be attributed to this effect, but, looking 
at the shape of the spectra in Figures 4 and 11, and as said be- 
fore, nonlinear mode coupling to the fastest growing modes, 
which here turn out to be at kjk\ k> 25 and 70, respectively, 
is a likely explanation. 

In this nonlinear stage, the ratio of parallel to perpendicu- 
lar energy begins to grow, and the field attains larger scale and 
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becomes more isotropic, as can be seen in Figure 1 . The a 2 
dynamo can be interpreted as an inverse cascade, in which 
parallel and anti-parallel magnetic flux is generated by the 
stretching of perpendicular flux, while the anti-parallel flux 
is kinematically "pumped" out of any given finite region of 
size L at a velocity of order a by the a effect. The pump- 
ing of flux into large regions of size L can be viewed from a 
modal point of view as inverse cascading of energy into small 
wavenumbers of order n/L. 

The observed scaling of the time evolution with J is rea- 
sonable: the timescale for J = 800 is about a factor of 
5.5 — lO 75 times less than that for J = 80. Because the 
stirring force of the counter-CR current is proportional to J, 
one might expect velocities to scale as J for a fixed correla- 
tion time. However, if an acceleration a operates over a set 
distance s, then the velocity scales only as (sa) 1 / 2 . Once 
the turbulent velocity exceeds vaq, the amplitude of the trans- 
verse motion associated with a given mode is <~ 2ir/k, and 
the correlation length varies more weakly than the correlation 
time with J. So we expect that the timescale for the dynamo 
mechanism scales as J s , where 0.5 < 5 < 1, and this is what 
we observe; S ~ 0.75. We have not verified how this scaling 
law extends to larger J. 

In Kolmogorov turbulence, the turbulent kinetic energy 
/ Ex(k) dk at any instant is proportional to P 2 / 3 , where P 
is the stirring power. In a situation such as the present one, 
where the stirring force F is proportional to JB, for a given 
stirring scale k, the velocity U scales as (F/k) 1 ' 2 , the power 
is proportional to F 3 / 2 and the energy should then, by di- 
mensional arguments, have a finite, steady state value that is 
proportional to F. The point is that for a given B, the energy 
E has a finite value to which it should rise and flatten out. 
There is indication of this in our simulation results, as seen in 
Figs. 2 and 3, where the total turbulent energy flattens out at 
toAofci = 0.45 for J = 80 (tvAoh = 0.08 for J = 800). 
It does not completely flatten out though, and we attribute 
this to the fact that £? rms is still creeping up with time due to 
the dynamo effect. Our simulations are not long enough to 
determine whether the magnetic energy Em always reaches 
equipartition with the kinetic energy Ek- If it does, then, be- 

1/2 

cause B scales as E M , dimensional arguments suggest that 

1/2 

the force F, which scales as JB, therefore scales as JE^ . 
So, if U scales as (F/k) 1 / 2 and thus Ek scales as F, then the 
force scales as J 2 . We believe it would take longer simula- 
tion runs that were presently feasible to check this prediction, 
and whether Em indeed scales as Ek- 

In the above discussion we have invoked the Bell instabil- 
ity to generate magnetic flux that is perpendicular to the local 
background, because this, from the point of view of the fluid, 
is an MHD effect and can be represented in an MHD sim- 
ulation. We note, on the other hand, that non-MHD effects 
could achieve the same result. For example, the non-resonant 
firehose instability could achieve the same stretching. The 
helicity that is required for the a effect relies on preferential 
growth of one circular polarization over the other, so the fire- 
hose instability on large scales, which has no such preference, 
could not by itself bring about an a effect. However, it could 
combine with the resonant cyclotron instabilities to do so. In 
this case, the firehose instability would play the role that dif- 
ferential rotation plays in the ail dynamo - that of creating 
a perpendicular field component, while the helical turbulence 
that is generated by the Bell instability would play the role of 



helical turbulence that, in the aft dynamo, is generated by the 
combination of convection in a stratified medium and Coriolis 
force. 

5. ASTROPHYSICAL APPLICATIONS 

In a real astrophysical system, there is a limited amount of 
time available. In the case of an expanding blast wave, this is 
of order the expansion time. Similarly, in an accretion shock, 
accreting matter continuously sweeps magnetic field down- 
ward over the crossing timescale of the accreting matter. The 
question is whether this is enough for significant magnetic 
field amplification. 

5.1. Blast waves 

In our simulations, the growth up to maximum amplitude 
takes on the order of 0.45 (0.08) Alfven crossing times across 
the box for J = 80 (J = 800); see Figs. 1 1-14. How- 
ever, this total time interval witnesses a gain of many orders 
of magnitude of the level of large-scale magnetic field, be- 
cause the simulation begins at the noise level. For J = 80 
it can be seen from Figs. 2 and 4 that at v^kit <~ 0.3, the 
growth at large scales speeds up, the energy in field pertur- 
bations at any given scale grows by a factor of <~ 10 3 / 2 per 
At <~ 0.04«Aofci, and, at somewhat later times, by a factor 
of nearly 10 2 per At <~ 0.04tA.ofci, suggesting that, during a 
period of exponential growth, the gain factor over an interval 
T at wavenumber k, is G(T, k) = io 2 ^ fe i' I 7 - 04 . For a su _ 
pernova remnant of radius R, the precursor has a width W of 
somewhat less than R, and the largest mode has a wavenum- 
ber of order ki = 2n/W > 1. Thus, the available time T for 
field amplification by cosmic ray current is T — W/u s ~ 
27r/(/cii>A-Myt)> wnere Ma = u s /va is the Alfven Mach 
number of the blast wave. Then vpJt\T <~ 2tt/Ma, suggest- 
ing that the a dynamo effect could amplify the field energy of 
the remnant by a gain factor G of order G = I0 7r / 01M ^ , and 
the amplification of the field's magnitude would be the square 
root of this factor, i.e. a factor of several for Ma ~ 10 2 . 
This appears to be consistent with Figure 13. According to 
Equation (1), and assuming P cr oc u 2 , we see that raising the 
potential CR current by ( 3 , speeds up the time evolution of 
the field amplification by ( s , where S is between 0.5 and 1; 
see the discussion above. So, very young SNR, where J can 
be much higher than in the runs we made, the gain factor for 
an expansion time could be higher. In Figure 4, for example, 
the magnetic energy on any given scale, [EK{k)k}/(k\/k) 
reaches, but does not significantly exceed, unity. 

For young, expanding supernova remnants Ma is typically 
10 2 , so the growth factor on the scale of the supernova could 
be of order one e-fold per expansion time if the assumptions 
of the simulation were to remain valid over the full interval. 
While not dramatic, neither is it insignificant in view of the 
uncertainties, a modest amount of magnetic field amplifica- 
tion may take place within an expanding supernova remnant, 
and this may be enough to be compatible with observational 
inferences (Pohl et al. 2005). While the correlation length of 
the magnetic field increases with time, it is at all times less 
than the size of the box by a factor of several. 

On the other hand, increasing the shock velocity from 
10~ 2 c to near c would raise J by a factor of 10 6 , thereby 
speeding up the rate of field amplification by J 6S > 10 3 
while decreasing the expansion time by only a factor of 10 2 , 
so there would then seem to be enough time to amplify the 
magnetic field by many orders of magnitude even on the 
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largest scale - as demonstrated in our simulations - if there 
are no other fundamental limitations that we have not yet 
identified. Thus, GRBs, which create ultrarelativistic shocks 
in the interstellar medium, would have the greatest potential 
for magnetic field amplification, because their ratio of particle 
pressure to initial magnetic pressure is so large. This is now 
discussed below. 

In our simulations, the velocity of the turbulence attains a 
magnitude of about 10 2 «ao, the magnetic field increases to 
~ 10 times the original field B , and the correlation length L 
of this field is about 10 _1 of the box size R. Thus the quan- 
tity eBL is not much changed from the original value eB^R, 
meaning that the maximum energy attainable by cosmic rays, 
<~ u s eBL/c (Eichler & Pohl 201 1), is not much changed by 
the magnetic field amplification. (However, until we are cer- 
tain how the final correlation length scales with the running 
time of the simulation, this matter remains not completely set- 
tled.) 

If the energy of the highest energy cosmic rays E is taken 
to be u s B R/c then the deflection of these highest energy 
implied by (though ignored in the simulations) is probably 
small for u rms «ii s . In general, the turbulent rms velocity 
u rms should be somewhat less than the shock velocity, u s , so 
the potential drop eu Yms BL/c across one correlation length 
L is somewhat less than the maximum energy attainable with 
shock acceleration eu s BL/c. This means that deflection is 
not a problem for the highest energy cosmic rays, but would 
be a problem for lower energy CR, and, to compute the cos- 
mic ray current, we are entitled to figure in only CRs of the 
highest energies, i.e. above u lms E max /u s . Magnetic field am- 
plification that uses CRs at lower energies is unlikely to help 
increase the maximum energy to which CRs can be acceler- 
ated by shocks. 

As mentioned above, the correlation length L is probably 
considerably less than the size of the box, i.e. the radius of the 
supernova remnant when applied to that context. It is possible 
that the rather large values of magnetic fields that have been 
claimed for young supernova remnants (10 3 G, i.e. about 10 2 
times the interstellar field of the Galaxy) are observationally 
compatible with such a small scale. 

5.2. Weakly Magnetized Relativistic Shock Waves 

We define a weakly magnetized relativistic shock wave as 
one where the kinetic energy of the upstream fluid flowing 
into the shock greatly exceeds the magnetic energy, i.e. where 
the upstream fluid is sufficiently magnetized as to be describ- 
able by MHD. Particles reflecting off the shock have velocity 
(3 S in the frame of the shock, where u s = (3 s c is the shock ve- 
locity which we assume to be equal to the streaming velocity 
discussed in Section 1. Assuming (1 — (3 S ) <C 1, velocities of 
up to [3 S + (1 — (3 s )/2 in the lab frame. The thickness of the 
shock precursor in the lab frame is thus of order i? s (l — /3 s )/2, 
where R s is the radius of the shock. The current density in re- 
flected ions, which we assume extend further upstream than 
the reflected electrons, is thus of order 

J~2en cT 2 s . (42) 

For typical GRB parameters n ~ 1cm -3 , u s ~ c, r s 2 = 
T s /100 (where T s is the shock Lorentz factor), Bo ~ 3/xG, 
the dimensionless parameter defined in the introduction is 
J ~ n Q m t c 2 /(B 2 /87:) = r 2 c 2 /v A 2 ~ 10 13 T 2 2 for scales 
fc _1 of order the reflected ion gyroradius. This is a far larger 
value than anything that can be reliably simulated, because the 



fastest growing mode is of too small a scale to be resolvable 
numerically. 

The thickness of the precursor of reflected ions from the 
shock is determined by how far ahead of the shock the ions 
can get before they are overtaken by the shock. A reflected 
ion by definition moves faster along the shock normal than 
the shock itself at the moment the ion crosses upstream, and 
it is overtaken by the shock after it has gyrated approximately 
l/r s of its gyroradius r g ~ T 2 miC 2 / 'eB , at which point its 
motion along the shock normal is less than that of the shock, 
so the shock overtakes it. Thus, it has moved a distance of 
Ar = r g /T s = r s m,iC 2 /eBo, whereas the shock has moved 
by f3 s Ar. The thickness of the precursor is then (1 — (3 S ) Ar ~ 
miC 2 /r s eB Q . 

The time over which a parcel of fluid at radius R is exposed 
to the reflected ion flux is 

At = Ar(l - p s )/c, (43) 

and its ratio to the Alfven crossing time across the precursor 
is 

= v A /c. (44) 

R/v A 

Due to numerical limitations, it remains unclear how much 
magnetic field amplification can take place. While we have 
shown that perturbations can grow by many orders of mag- 
nitude, we do not have any runs in which the final magnetic 
field was more than a factor of 10 or so more than the original 
field. As we do not understand the nonlinear dissipation of 
magnetic field, so we do not know a priori how the limiting 
field scales with J. Yet we argue on theoretical grounds that 
the a-dynamo should be able to amplify the field at relativis- 
tic or near relativistic shocks by a large factor, as it apparently 
does in many compact objects, and the existence of such a 
dynamo is ultimately due to the left-right asymmetry of the 
magnetic turbulence that is generated by the cosmic rays. 

5.3. Cosmic Rays in the Galaxy 

Let us now consider systems with lifetimes that are large 
compared to the Alfven crossing time, such as the Galaxy. 
Here there is enough time for the standard ail mean-field dy- 
namo dynamo to work if there is a source of right-left asym- 
metric turbulence (Ruzmaikin et al. 1988; Brandenburg & 
Subramanian 2005). In any system where the magnetic en- 
ergy has attained rough equipartition between magnetic field 
and cosmic ray pressure, the ratio on the right hand side of 
Equation (1) obeys 87rP cr /_B 2 < 1. The first ratio on the 
right hand side of Equation (1) u s /c, must also be less than 
unity. The third ratio in Equation (1) must also be less than 
unity for the assumption for the Bell instability (small de- 
flection of the current-bearing cosmic rays) to be valid. It 
follows that for the above assumptions, the left hand side 
J = 47r J CI /ckiBo -C 1. This, however, means that the Bell 
instability does not occur, as seen from Equation (9); rather, 
the effect of the cosmic rays is to slightly decrease the phase 
velocity of stable Alfven waves. It follows that Bell turbu- 
lence cannot amplify the magnetic field to near equipartition 
with the cosmic ray pressure. Some estimates in the literature 
use a cosmic ray density of <~ 10 9 cm~ 3 . It would follow 
that J > 1. However, the problem here is that this value 
for the cosmic ray density includes low energy (GeV) cosmic 
rays, which satisfy the assumption for small deflection only at 
extremely small spatial scales. It is doubtful that such modes 
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should even be Bell-unstable at all because of the ~r\k 2 damp- 
ing term in the imaginary part of the frequency, as expressed 
by Equation (9), which grows with k faster than the growth 
term, and which should therefore dominate at small spatial 
scales. 

As discussed in the introduction, the factor 3P cr / (B 2 /An) 
is of the order of unity for the Galaxy as a whole if all relativis- 
tic CRs are included. The anisotropy u s /c is bounded by ob- 
servations to be at most 10 ~ 3 . Finally, the term eB / kTrriiC 2 
must be less than unity to satisfy the small CR deflection cri- 
terion that is the basis for the Bell instability. Altogether it 
follows that the quantity J in Equation (1), is less than unity 
in the Galactic disk. This, however, implies Bell stability by 
Equation (9). So, whereas the standard a£l mean-field dy- 
namo (Ruzmaikin et al. 1988; Brandenburg & Subramanian 
2005) and tangling the Galactic magnetic field with cosmic 
ray streaming instabilities may be viable ways to amplify the 
Galactic magnetic field, resonant CR streaming instabilities 
seem the more promising way to do it over large volumes, 
where 3P cr / (B / 4tt) is of the order of unity. Resonant CR 
streaming instabilities produce Alfven modes of preferential 
circular polarization, just as the Bell instability does, so the 
theoretical mechanisms for field amplification that are dis- 
cussed in this paper apply to them as well. 

We conclude that dynamo activity in the Galactic disk from 
the twisting of the field by collective cosmic ray interactions 
can take place in principle. The numbers seem marginal, so 
more careful investigation is needed to settle this point. 

5.4. Unconfined intergalactic cosmic rays 

In intergalactic space, the cosmic ray pressure is compara- 
ble to the magnetic pressure and there is then the chance that 
J greatly exceeds unity. Assuming that streaming instabili- 
ties keep the streaming velocity below the Alfven velocity va, 
we can replace u s with va, so in principle the field can be am- 
plified to the point where the magnetic pressure B 2 /4n is of 
the order of P ci va / c. 

We use the following parameters for plasma and cosmic ray 
particles: we assume u CI = 3 x (10 6 -10 7 ) cms -1 for the drift 
velocity of cosmic ray particles, n CI = 10~ 9 -10~ 10 cm~ 3 
for the number density of cosmic rays particles, rn = 10~ 4 - 
10~ 2 cm~ 3 for the mean number density of plasma ions, 
B* = 1 fiG for the equilibrium mean magnetic field, and 
£o = 3 x 10 18 cm for the turbulent length scale. The Alfven 
speed is va — 10 8 -10 9 cms -1 and the dimensionless J pa- 
rameter is 



which is comparable to the values studied in this paper. 

6. CONCLUSIONS 

We have investigated the mean-field dynamo mechanisms 
in a turbulent plasma with a cosmic ray current. We find 
a linear growth stage, corresponding to the Bell instability, 
and then a nonlinear stage corresponding to the production of 
a fully-developed MHD turbulence and generation of larger- 
scale magnetic field. 

In the nonlinear stage, the level of MHD turbulence con- 
tinues to grow (see Figs. 2 and 4), and the correlation length 
increases with time (see Figure 9). The turbulence develops 
a Eu{k) oc fc~ 5 / 3 spectrum, despite the fact that the initial 
spectrum is EM{k) oc fc 4 . These effects were far from obvi- 



ous, even given the stirring by the counter- J cr x B forces on 
the fluid, which is the basis of the Bell instability. 

We suggest that this combination of magnetic field ampli- 
fication and large-scale ordering is due to the a 2 dynamo 
instability, in which the original field is stretched by unsta- 
ble Alfven modes, and then the stretched component is it- 
self stretched by unstable circularly polarized Alfven modes. 
This effect is to be contrasted with the long wavelength linear 
instability discussed by Bykov et al. (2011a), because their 
growth rate vanishes when the perturbations are independent 
of z. According to the simulations, the level of magnetic field 
on any scale can be enhanced by a factor of several within 
one expansion time, even at the largest scales, but there is no 
direct numerical evidence at present that it can be enhanced 
by much more than that. Nevertheless, the analysis presented 
here predicts that much larger enhancement, via exponential 
growth to the a effect, is possible. 

We have also used DNS and the test-field method to con- 
firm our analysis. DNS shows that the instability has three 
stages. In the first stage, the Bell instability is excited; in the 
intermediate stage the linear growth continues among the high 
k modes, while mode coupling feeds the low fc, large-scale 
modes. In the third stage, growth on large scale continues, 
apparently due to the a 2 dynamo. 

We find, as expected, that the value of a, which has units 
of velocity, can never be much greater than the rms turbulent 
velocity, u rms . The rms turbulent velocity, u rms , by energy 
conservation, must be less than the shock velocity u s , prob- 
ably much less. Thus the maximum scale to which the a 2 
dynamo can operate efficiently must be limited to a/u s <C 1 
of the radius of the blast wave R; because a < u Tms <C u s , 
there is not enough time in a single expansion time to stretch 
or move the field significantly over the characteristic expan- 
sion time R/u s of the blast wave. There is enough time only 
to greatly amplify the field on scales much smaller than R. 
Similar remarks would apply to accretion shocks, since ma- 
terial is swept out of the region of field amplification over a 
timescale of R/u s . 

In blast waves from GRBs, on the other hand, this maxi- 
mum scale of amplification, aA/c (where A is the scale of 
the cosmic ray precursor), could be much larger than the ion 
skin depth, so dynamo activity in the precursors of such blast 
waves could significantly increase the magnetic correlation 
length relative to what is produced by the Weibel instability, 
so that the Ohmic dissipation downstream would not be so 
devastating. 

The maximum energy to which cosmic rays can be acceler- 
ated by expanding supernova remnants, which is proportional 
to B R without magnetic field amplification, is not greatly 
enhanced by large-scale field amplification (not, at least, in 
the parameter regime we have explored), since the increase 
in field strength comes at the cost of decreased correlation 
length L, and the product BL is enhanced far less than B. 
However, quantification of this point, made in greater detail 
by Eichler & Pohl (201 1), would require a more extensive li- 
brary of simulations, for all values of the ion current parame- 
ter J that could conceivably occur in nature, and such simula- 
tions become very difficult for high J. Clearly if the pressure 
of initial field is sufficiently small relative to the cosmic ray 
pressure, then extremely large values of J are possible, and 
the range J ^> 10 3 has not been explored numerically. On 
the basis of extrapolation of the two values of J (J = 80 
or 800), for which we run high-resolution simulations, the e- 
folding timescale for field amplification t(J, fc) over scale 
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Trfc- 1 is approximately C^AfcJ 5 ] -1 ; (C ~ 1,5 ~ 0.7), 
while we expect on the basis of Equation (1) that the largest 
allowable value of J, which occurs at eB /[kTimc 2 ] <~ 1, 



scales as 



as v 



-1+25^-1 



So the timescale for field amplification scales 
Clearly this timescale 



can be made small enough for sufficiently weak fields and 
sufficiently small spatial scales (i.e. high k). How field am- 
plification at these scales and field strengths would ultimately 
affect the capability of shock acceleration remains an impor- 
tant question for future research. 

We have also considered the potential of this dynamo pro- 
cess in steady situations where newly created magnetic flux 
does not quickly sweep through the region where cosmic rays 
stream. An example of this could be the Galaxy itself, with 
escaping CR providing a steady flux over the lifetime of the 
Galaxy. If we assume that the CRs that provide current for 
growth on scale fc _1 must satisfy eB /[kTrriiC 2 } < 1, then, 
by Equation (1), the condition that J > 1 requires that 
P CI » (B 2 /8tt)(c/u s ) (where u s now stands forCR stream- 
ing velocity), so the dynamo can only bring the magnetic pres- 
sure to within a fraction u s /c of the CR pressure. On the other 



hand, this may be interesting, because it could "sprout" a seed 
magnetic field to a sufficiently large amplitude that some other 
mechanism, such as the magnetorotational instability, could 
further raise the field pressure to its present value, <~ P CI . 

We have also noted that resonant CR streaming instability, 
coupled with differential rotation, remains a possible way to 
promote an aSl dynamo in the Galactic disk. 
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APPENDIX 

A. DERIVATION OF EQUATION FOR TOTAL a EFFECT 



We determine the contributions to the mean electromotive force, £(B) = uxb, caused by cosmic ray particles in homogeneous 
turbulent plasma. The procedure of the derivation of equation for the mean electromotive force is as follows. We use a mean field 
approach in which the magnetic and velocity fields are divided into the mean and fluctuating parts, where the fluctuating parts 
have zero mean values. The momentum and induction equations for the turbulent fields are given by 



. du(t, x) 



■Vptot + -!-[(&• V)B + (B-V)b] - -J™xb+-en CI (uxB) + N u , 

47T C C 



db{t dt X) = (B-V)u - (u-V)B - B (V-u) + N b , 



(Al) 
(A2) 



where u and b are fluctuations of velocity and magnetic field, B is the mean magnetic field, J cr is the mean density of the 
electric current of cosmic ray particles, N u and N b are the nonlinear terms which include the molecular dissipative terms, 
Ptot — P + (B-b)/4:ir are the fluctuations of total pressure, p are the fluctuations of fluid pressure. To exclude the pressure term 
from the equation of motion (Al) we calculate Vx(Vxtt). Then we rewrite the obtained equation and Equation (A2) in Fourier 
space. 

A.l. Two-scale approach 
We apply the two-scale approach, e.g., a correlation function, 



•r + iK-R) 



m(x)uj(y) = J dkidk,2 Ui{k\)uj{k2) exp{i(fci-x + &2 - y)} = J dk dK fij(k, K) cxp(ik- 
= J dk fij(k, R) exp(ife-r), 

(see, e.g., Roberts & Soward 1975). Hereafter we omitted argument t in the correlation functions, (fe, R) — L(m,; uj), where 



L(a; c) 



a(k + K/2)c(-k + K/2) cxp (iK-R) dK, 



and we introduced new variables R = (x+y)/2, r = x—y, K = fci+fe2,fc = (fei— fe2)/2. The variables R and K correspond 
to the large scales, while r and fc correspond to the small scales. This implies that we assumed that there exists a separation of 
scales, i.e., the maximum scale of turbulent motions £o is much smaller than the characteristic scale Lb of inhomogeneity of the 
mean magnetic field. 

A.2. Equations for the second moments 

We derive equations for the following correlation functions: fij(k,R) = L(in;Uj), hij(k,R) = (4:ir~p)~ 1 L(bi] bj) and 
gij(k 7 R) = L{bi] Uj ) . The equations for these correlation functions are given by 

dfn(k) 



dt 



= i{k-B)$ij + D im (ki)f mj + D jm (k 2 )f im + A 4m (fci)5 mi (fc) + A im (k 2 )g im (-k) + //• + ffi, 



(A3) 
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= -i(fc-B)$ y + ifc„bi„(fe) - flj -„(fc) + J* + (A4) 

= i(k.B)[f l3 (k) - h l3 (k) - fc<f >] - ik m B t f mj + D jm (k 2 )g im (k) + (4ttp) A jm (fc 2 )/ im + J» + (A5) 

where <J>y (fc) = (47rp) _1 [^ (fe) — gji(— k)]. Hereafter we omitted argument R in the correlation functions and neglected terms 
~ 0(V 2 R ), 



n. . _ o P . . ncrjL 4..— o P .. 7 cr k O cr — 7 cr _ fAfi"> 



_ j cr = 

2c/5 ' 2c/?' 



£y„ is the fully antisymmetric Levi-Civita tensor, the terms fW, hfj and gf^ are determined by the third moments appearing due 
to the nonlinear terms, the source terms l{- , 1^ and I?- which contain the large-scale spatial derivatives of the mean magnetic 
and velocity fields, are given by Eqs. (A3)-(A6) in Rogachevskii & Kleeorin (2004). These terms determine turbulent magnetic 
diffusion and effects of nonuniform mean velocity on mean electromotive force. In the present study we neglect small effects of 
cosmic ray particles on the turbulent magnetic diffusion. 

For the derivation of Eqs. (A3)-(A5) we use an approach that is similar to that applied in Rogachevskii & Kleeorin (2004). We 
took into account that the terms with symmetric tensors with respect to the indexes i and j in Equation (A5) do not contribute 

to the mean electromotive force because £ m — e m ji gij. We split all tensors into nonhelical, hij, and helical, h[f \ parts. The 

(H) (H) 

helical part of the tensor of magnetic fluctuations h\ - ' depends on the magnetic helicity, and the equation for h\ - ' follows from 
the magnetic helicity conservation arguments (see, e.g., Rogachevskii & Kleeorin 2004; Brandenburg & Subramanian 2005, and 
references therein). 

A.3. r-approach 

The second-moment equations (A3)-(A5) include the first-order spatial differential operators AT applied to the third-order 
moments M<- m K A problem arises how to close the system, i.e., how to express the set of the third-order terms AfM^ lir> through 
the lower moments M^ u \ We use the spectral r approximation which postulates that the deviations of the third-moment terms, 
AfM( IU ) (fe), from the contributions to these terms afforded by the background turbulence, MM^ 111 ' ^ (k), are expressed through 
the similar deviations of the second moments: 



A^M( m )(fe) -^"'"'(fc) = j- [MW(fe) - M( n '°)(fc) 

T(k) L 



(A7) 



(Orszag 1970; Pouquet et al. 1976; Kleeorin et al. 1990; Rogachevskii & Kleeorin 2004), where r(k) is the scale-dependent 
relaxation time, which can be identified with the correlation time of the turbulent velocity field. The quantities with the superscript 
(0) correspond to the background turbulence (see below). We apply the spectral r approximation only for the nonhelical part hij 
of the tensor of magnetic fluctuations. 

A.4. Solution of equations for the second moments 

First we solve Eqs. (A3)-(A5) neglecting the sources l(- , 1^ , with the large-scale spatial derivatives. Then we will take into 
account the terms with the large-scale spatial derivatives by perturbations. We subtract from Eqs. (A3)-(A5) the corresponding 
equations written for the background turbulence, use the spectral r approximation and neglect the terms with the large-scale 
spatial derivatives. We assume that the cosmic ray velocity is much larger than fluid velocity, so that the terms oc in 
Eqs. (A3)-(A5) vanish. Next, we neglect the effect related to the compressibility of the turbulent velocity field. Such effects are 
important when the Mach number is of the order of or larger than 1 . We also assume that the characteristic time of variation 
of the second moments is substantially larger than the correlation time r(k) for all turbulence scales. This allows us to get a 
stationary solution for the equations for the second-order moments, M^ ll \ Thus, we arrive to the following steady-state solution 
of Eqs. (A3)-(A5): 



(A8) 



fij(k) m f\f{k) + ir(fc-B)$«(fc) + r A im g m j(k) + A jm g im (-k) 

hij(k) w hfik) - ir(k-B)^ij(k), (A9) 

9ij (k) w glfik) + ir{k-B) [,4(fe) - fcy(fc)] + t (47rp) A jm h im (k), (A10) 

where /y , hij and g^ are solutions without the sources l(- , 1^- and Ifj . In the present study we consider linear effects in 
perturbations of the mean magnetic field. The nonlinear mean-field modeling in turbulent compressible MHD flows with cosmic 
rays is a subject of a separate ongoing study. 

A.5. Model for the background turbulence 
Now we need a model for the background anisotropic turbulence [see Eqs. (A8)-(A10)]. The anisotropy is caused by the 
equilibrium mean magnetic field B„. Generally, a model of an anisotropic turbulence with one preferential direction can be con- 
structed as a combination of three-dimensional isotropic turbulence and two-dimensional turbulence in the plane perpendicular to 
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the preferential direction (see, e.g., Elperin et al. 2002). Also we take into account that the tensor (fc) is the sum of non-helical 
and helical parts of the turbulence. A non-zero kinetic helicity is caused by the Bell instability. To relate the velocity fluctuations 
tensor f^\k) with the magnetic fluctuations tensor hf^(k) and the cross-helicity tensor g^\k), we use the relation between 

the magnetic and the velocity fields in the Bell mode: f/°)(fe) = i(k-B*) u^°\k)/^B, where 75 is determined by Equation (9). 
We use the following model for the background anisotropic homogeneous and helical turbulence caused by the Bell instability: 

/f(fc) = ui 0) (*) U f(-fc) = g| I [(1 - e) (5« - hj) + 26 (S tJ e iej - k±)] WW - ^ e ijn k n ««>) . (V x ««»)} , 
(All) 

,(o )(fc) s ^j^HO = (Lcr fc) ; (o )(fe); (A12) 

^(fc)^^(fc) U f(-fc) = (4 7 rp^fc) 1/2 *^fS\k), (A13) 

where L cr = c£?*/(47rJ cr ), £?(fc) = (g — 1) £0 (^0 k)~ q is the energy spectrum function, the length £q is the maximum scale of 
turbulent motions, e is the unit vector directed along the equilibrium mean magnetic field £?„, Sij is the Kronecker unit tensor, 
< e < 1 is the anisotropy parameter of turbulence, k = fc + k z e, k z — (k-e), kij = kikj/k 2 and k^j = kj-kj- / (fc" 1 ") 2 . The 



turbulent correlation time is r(fc) = C T r (£0 k) ^, where the time r = £o/u , u = y [it(°)] 2 is the characteristic turbulent 
velocity in the scale £ n and C T is the coefficient. For the background turbulence with a constant dissipation rate of turbulent 
energy in inertial range of scales, the exponent fi = q — 1, the energy spectrum E(k) oc —dr/dk and the coefficient C T = 2. 

Using the solution of the derived second-moment equations (A8)-(A10), we determine the contributions to the mean electro- 
motive force, £f = £i mn J b n (k) u m (—k) dfc, caused by cosmic ray particles in homogeneous turbulent plasma. 

A.6. Derivations of contributions to the a effect 

We take into account effects which are linear in the perturbations of the mean magnetic field: B = B — £?*, i.e., we consider 
a kinematic mean-field dynamo. Substituting Eqs. (A8)-(A9) into Equation (A10) we obtain: 

9ij (k) + £ I] (k) + £ H \k), (A14) 



g^(k)^ir(k-B) [fi;'(k)-h^(k) 

(k) « zr 2 (k-B) [A rm (fc) + A jm g<& (-fe) - (Anp) A jm S<°> (fc) 
= ir 2 {k.B)\A im g^){k) 3A jm g^(k) 



mru\ _ £(0) 



(A15) 



(A16) 



gS"\^g§\k) + r (47rp) A jm h$\k), (A17) 

where we have taken into account that g^\k) = g^\k) = — #-^(— fc). The mean electromotive force is given by £" — 

£imn J g n m (fc) dfe, where the tensor cnj (fc) is determined by Equation (A14). There are two contributions to the a effect caused 
by: 

(i) non-zero kinetic helicity produced by Bell-instability; this contribution is determined by the tensor g^\k), where the 
background turbulence (fc) is described by the term, oc — (i /fc 2 ) £ij n k n u^> ■ ( V x ) in Equation (Al 1); and 

(ii) interaction of the mean cosmic ray current with small-scale anisotropic turbulence; this contribution is determined by the 
tensor g^j T \k), where the background turbulence f-j\k) is determined by the term, oc e (5^ — e,ej — fc^) in Equa- 
tion (All). 

The first contribution to the mean electromotive force caused by a non-zero kinetic helicity effect in the Bell turbulence is given 
by 

£P = e mn J gW(k) dfc = ie imn J r(fc) (k-B) (fc) dfc = a$ Bj, (A18) 

where a\p — a™ 5^ and af is determined by Equation (27). In the derivation of Equation (A18) we have taken into account 
that L cr ki = J^ 1 <C 1 that allows us to drop contributions oc hff (fc) in comparison with that proportional to f-^ (fc). 
The second contribution to the mean electromotive force, £ • !I \ is caused by the non-helical part of the anisotropic turbulence: 



£ (II) -e 



, J (fc) dfc - 4 i e mrn J t 2 (fc) (k-B) A mp g<$ (fc) dfc - a<f> B 3 , (A19) 
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where the tensor — a" (Sij+eiCj) and a" is determined by Equation (14). The third contribution to the mean electromotive 
force, £j- 7/7 ^ = Si mn J gnm\k) dfe, is constant and, therefore, does not affect the large-scale dynamo. 

A.7. Integrals used in Section A.6 
To integrate over the angles in fc space we used the following identities: 



J kij n sgn(fc 2 ) sin 9 d9 dip 



, 3 — sgn(fc z ) sin 9 d6 dip = — 
k ' 2 



(A20) 
(A21) 



where Pi„(e) = <5y — e^ej and k z = k cos 6*. 



A. 8. 77?e realizability condition 

Let us consider the case when the spectral functions for the kinetic helicity, x{k), and turbulent kinetic energy, Uq E(k), are 

different, where ix(°) • (V x «( )) = J \(k) dk and [it(°)] 2 = u^J E(k) dk. The realizability condition for the kinetic helicity 
(Moffatt 1978) reads: 

X (k) < ul E{k)k. (A22) 

Let us determine the explicit expression for the function x(k) using Equation (26) for the estimate for the the kinetic helicity 
tt(o). x u,(°)) for the Bell background turbulence: 



«(°)- (V x 



T JJ 



4?r J c 



1/2 



2 r 

Va^ J (e kf /2 T(k)E(k)dk, (A23) 



where = (q — 1) (^o fc) _<? an d r (^) = 2 To (i'o A;) 1-9 . Therefore, the function x(fc) is given by 

x(fc) = Vjv A ^(e kf 2 - q E(k), 

and the realizability condition for the kinetic helicity yields: 



u 



where 



4tt J cr £ 



(A24) 



(A25) 



(A26) 



and the Kolmogorov spectrum corresponds to q = 5/3. 
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